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1. Introduction 

One of the main aspects to be studied in thermal quantum chromodynamics (QCD) is the 
transition - or smooth crossover - from a chirally symmetric and deconfined region of phase space 
to a hadronic phase with broken chiral symmetry. This transition has occurred in the cooling of the 
early universe and is hoped to be reproduced in heavy ion collision experiments. 

A lot of effort has been put into lattice studies of this transition, see e. g. [|l|] and the most recent 
overview by K. Kanaya [Q]. Hybrid Monte-Carlo simulations are restricted to vanishing or at most 
small chemical potential due to the sign problem that otherwise prohibits importance sampling. 
In case of zero chemical potential and physical quark masses, it has been demonstrated that the 
thermal transition is no true phase transition but really a smooth crossover [Bp. The current knowl- 
edge about the nature of the phase transition for different quark masses is shown in Figure 0. This 
dependence is particularly important for lattice simulations with Wilson-type fermions, for which 
physical quark masses have not been feasible so far and extrapolations are necessary. Furthermore, 
by knowing the qualitative properties of this plot, some constraints on the enlarged phase diagram 
including a non-vanishing chemical potential can be anticipated. 



The critical temperature for a true phase transition is uniquely determined by non-analyticities 
in the partition function. For a crossover there are only pseudo-critical temperatures which depend 
on the observable. Here we use the Polyakov loop and its susceptibility and compare to the pion 
norm and the screening mass in the pseudo-scalar channel. 

The aforementioned lattice results have been obtained from simulations in which predomi- 
nantly staggered fermions have been applied. However, the question to what extent this type of 
fermions is reliable has still not been completely settled. The main problem is the so-called rooting 
procedure, i. e. the need to take the fourth root of the staggered fermion determinant in order to 
describe a single flavour of quarks, see e.g. [Q]. Therefore, it is highly desirable to use different 
fermion formulations so that the possible caveats are under better mutual control. Wilson-type 
fermions have very different systematic effects from the staggered formulation. In particular, they 
do not suffer from the so-called taste-breaking that needs to be controlled in staggered simulations, 
see [j5j]. The major problems of Wilson-type fermions are of course the explicit chiral symmetry 
breaking and the introduction of ff(a) effects. 




Figure 1: Nature of the thermal transition in the quark mass plane. 
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In this contribution we report on our ongoing investigation of the thermal transition of QCD 
applying so-called maximally twisted-mass fermions. These Wilson-type fermions offer automatic 
@{a) improvement. The quark mass is determined by the multiplicatively renormalisable twisted 
mass parameter jX. For an introduction, see [||]. 

We use one doublet of degenerate fermions for the twisted mass action. Thus, in terms of 
Fig. [j], we are restricted to the upper boundary where the strange quark is infinitely heavy. Here 
the situation may well be different from that at the physical point. However, the Nf = 2 setup is 
somewhat natural to start from with our action as the twisting always acts on fermion doublets. 
Furthermore, there remains the question whether the thermal transition in the two-flavour chiral 
limit is really of second order in the 0(4) universality class. Therefore, the chosen setup is very 
well suited to establish further information about the quark mass dependence of the QCD transition. 

Besides our approach relying on the twisted mass formulation, also clover-improved Wilson 
fermions are currently used for finite temperature studies [0, §], see also [l(J in these proceed- 



ings. As compared to our previous work [Mil, 12], we now report on simulations with light quark 



masses > 300 MeV that will give information on the quark mass and lattice spacing dependence of 
the critical temperature. 

In the following section, we present the maximally twisted mass lattice formulation as used 
for our simulations. Furthermore, we describe the non-trivial (k",j8,/i) phase diagram and explain 
our realisation of scans in /3. Section ^] gives a summary of our results obtained so far. We show 
the Polyakov loop and the flavour non-singlet pseudo-scalar screening mass as examples for the 
observed signals and collect the transition temperatures available from our data. This also allows 
for a discussion of cutoff effects. Finally, in Section || we conclude and describe our plans for the 
continuation of this work. 

2. Lattice action 

For our simulations we rely on zero temperature information by the European Twisted Mass 
Collaboration (ETMC) [|l^, [Wj]. The gauge action is chosen to be tree-level Symanzik improved. 
The maximally twisted fermion matrix takes the form 

M mtm = 1 - K c D wi i son + 2i K c an y 5 T 3 (2.1) 

and acts on one fermion doublet. The Pauli matrix T 3 acts in flavour space and K c is the hopping 
parameter at its critical value. For jj, = the action has its standard Wilson form. In the continuum, 
the new term ia/i/sT 3 can be obtained by a non- anomalous chiral-flavour rotation. However, this 
rotational symmetry is broken by the Wilson-term so that one has different cutoff effects for each 
value of this so-called twist-angle. For maximal twist, i. e. if k = Kc, it has been shown that non- 
vanishing physical observables are automatically G[a) improved (see ||6|] and references therein). 



We have checked that this remains true at finite temperature by means of a quenched analysis [ |1 1| ] 
and a perturbative calculation []i"5|]. 

Wilson-type fermions suffer from unphysical phases in their bare parameter phase space. For 
standard Wilson fermions the Aoki phase at strong coupling is well known [[if], [17[ ]. For the en- 
larged (k,P,h) space of twisted mass fermions an additional phase transition plane has been found 
that incorporates the line of critical hopping parameters fQ-(/3,jU =0) t]T8[ |l9|,p 
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Figure 2: Left: Phase diagram in (k, j3,/i)-space for twisted mass lattice fermions at finite temperature. 
Right: Interpolation of JQ(j3) from ETMC results. 
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Figure 3: Real part of the Polyakov loop (left) and its susceptibility (right) from the two available lattice 
spacings at m K w 380 MeV. The vertical line indicates the value of T c . 



It is of course very important to make sure that the finite temperature simulations are not 
affected by these unphysical phases. In Fig. ||, left, our findings are summarised [12]. The thermal 
transition can be found as a surface that winds around the line of critical hopping parameters. This 
behaviour can be understood as a remnant of the continuum symmetry under twist rotations . 
In conclusion, it is possible to perform physical simulations towards the continuum limit that are 
not affected by unphysical bulk behaviour. 

Our simulation strategy is to perform scans in the lattice coupling j8 while keeping maximal 
twist in order to be G(a) improved. The knowledge of JC c (j8) as well as the lattice spacing and pion 
masses can be obtained from ETMC data [13, 14]. The interpolation of K" c (j3) is shown in Fig. ||, 
right. It is also possible to keep a fixed pion mass by adjusting ^i(P,m n ). 



3. Thermal signals 

Since the expected behaviour is that of a continuous crossover, it is very important to study 
the temperature dependence of several observables. In our study, we look at the real part of the 
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Figure 4: Left: Flavour non-singlet pseudo-scalar screening mass from the two available lattice spacings 
at m n ss 380 MeV. Right: Variance of pion norm at m% « 380 MeV. The vertical line marks T c from the 
Polyakov loop susceptibility. 



Polyakov loop, the chiral condensate, the pion norm (pseudo-scalar susceptibility), the plaquette 
and screening masses. 

In this contribution we present plots for the Polyakov loop, Fig. ||, the pseudo-scalar flavour 
non-singlet screening mass and the variance of the pion norm, Fig. |] at m % ss 380 MeV. The peak 
position of the susceptibility of the Polyakov loop is used to define T c . For m n 380 MeV, we 
obtain T c (N r = 10) = 243(10) MeV and T c (N r = 12) = 243(8) MeV from the respective maximal 
values. 

These crossover temperatures are consistent with the pseudo-scalar flavour non-singlet screen- 
ing mass, presented in Fig. [| left. This screening mass shows a qualitative change in behaviour at 
the onset of the transition region, i. e. just left to the corresponding value of T c . A combined anal- 
ysis of this screening mass and the corresponding scalar one will give information on the splitting 
introduced by the L^^-anomaly. If the anomaly is still sufficiently strong near T c , the transition 
in the chiral limit is expected to be in the 0(4) universality class, although a first order scenario 



is not excluded [22, 23]. The size of this splitting is thus an indicator for the order of the Nf = 2 
transition in the chiral limit. 

The variance of the pion norm normalised to the spatial volume 

of,p=^(((|»p) 2 )-(w 2 ) 2 ) (3.1) 

is shown in Fig. ^[ right. For this observable a peak can be found at approximately the same position 
at which the Polyakov susceptibility peaks. 

In Fig. ||, left, T c for the two lattice spacings at m n « 380 MeV is plotted. The N t = 12 result 
can possibly be further refined and at least a third lattice spacing is needed for a final statement. 
However, the cutoff dependence appears to be rather small. 

In addtion to the points for m n rj 380 MeV we currently simulate at N t = 12 with m n rj 
300 MeV and 460 MeV For m n « 300 MeV, a first value of T c w 225(11) MeV is available. Fig. |, 
right, shows our results compared to those of ^7§. It seems that our points tend to be at somewhat 
higher temperatures. A small shift in our results due to the scale set by the new ETMC data might 
occur, but we expect this to be at the percent level. 
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Figure 5: Left: Lattice spacing dependence for T c as detemined by the Polyakov loop for m n ps 380 Me V. 
Right: Values of T c (m K ) from the Polyakov loop susceptibility for different lattice spacings. The values of 
QCDSF/DIK from are shown for comparison. 



4. Conclusions and outlook 



We have presented results of our study of the thermal transition of QCD applying maximally 
twisted mass fermions. Relying on the ETMC data at zero temperature, we can investigate a mass 
range m % > 300 MeV. Our current simulations will soon provide enough input to perform a chi- 
ral and continuum extrapolation of the critical temperature. So far we have concentrated on the 
Polyakov loop to define this temperature. Screening masses supplement our set of observables. 
Their splittings can give insights into the Nf = 2 chiral transition. 
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